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Charge Density Wave in 
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Recent experimental results [ 1 , 2 1 revealed the stability 
of a peculiar charge density wave (CDW) with odd pe- 
riod (A = 3, 5) in the coupled CU2O5 spin ladders in 
Sri4_ x CaxCu2404i, while the CDW with even period is 
unstable. Here we resolve this contradiction by solving 
the effective t-J model for the two coupled chains sup- 
plemented by the intersite interaction term V. 



1 Introduction It is at the heart of the understanding 
of the high-temperature superconducting (SC) phase to in- 
fer what kind of state competes with the SC one in the pos- 
sible quantum phase transition between the two. As in the 
two-dimensional (2D) copper oxides one of the candidates 
of such a state is the stripe phase, it is natural to expect that 
in the quasi-2D copper oxide coupled-ladder system also a 
kind of charge ordered state would compete with the pre- 
dicted there SC state |4|. Indeed in Sri4_ a ;Ca a .Cu2404i, 
which has the SC state for x = 13.6 and under pressure 
larger than 3 GPa Q, a charge density wave (CDW) phase 
with period A = 3 and A = 5 was observed when the 
Calcium concentration was tuned to x = 11 (which cor- 
responds to the hole concentration rih = 6/5 per copper) 
and x = (where n/j = 4/3), respectively QE)- Although 
the charge transfer model for coupled ladders proposed in 
Ref. may explain the onset of such an odd-period CDW 
state, it is the t- J-like model which is a natural Hamil- 
tonian for the doped spin ladders COG]. However, then this 
experimental observation is a challenge for the theory since 
such a novel CDW state with the odd period has not been 
predicted by the standard t-J model for the ladders 
Moreover, the standard t-J model predicts the CDW order 
with period A = 4 to be stable for nu = 5/4 which would 
correspond to the x = 4 doping in Sr^-zCa^O^O^ 
where no charge order was observed 0]|2]. 



Crucial to this finding is the novel interaction term V 
(responsible for the interladder repulsion between near- 
est neighbor holes), which is derived from the on-site 
Hubbard-like repulsion between electrons on the same 
oxygen sites but within two distinct orbitals which be- 
long to two neighboring CU2O5 ladders in the charge 
transfer model 0. 
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The purpose of this paper is to shed some light on this 
apparent contradiction between the experiment and the the- 
ory. Therefore, first we discuss which terms are missing in 
the standard t-J model while in the second part of the pa- 
per we look at the generic features of this new model by 
solving it in a simplified case for two coupled chains using 
the Gutzwiller approach and the mean-field approximation. 

2 Model Our starting point is the effective one- 
dimensional (ID) t-J model for a single chain, 



) , (1) 



where: (i) the operator d\ a = d\ a (l — Tiid,-a) creates an 
electron at site i with spin a in the restricted Hilbert space 
with no double occupancies, (ii) hida = d\ a d ia is the cor- 
responding electron number operator, and (m) Si is an 
5 = 1/2 spin operator. Besides, the average number of 
d holes per site here is n = 2 — nu, with n = ~}2 a {nidcr)- 

In the standard approach to the CU2O5 systems this 
model is derived from the charge transfer model which is 
similar to that for CuC<2 planes (7J. Then the model should 
also contain the kinetic oc t and superexchange oc J terms 
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Figure 1 Artist's view of the three coupled Q12O5 lad- 
ders (denoted by grey, white, and again grey color, re- 
spectively). Only orbitals which are included in the charge 
transfer model of Ref. |3) are shown (where x = 2p x , 
y = 2p y , and d = id x 2_ y 2). The vertical box depicts two 
coupled CUO3 chains belonging to two different neighbor- 
ing ladders and discussed in the paper. 



along the rungs of the ladder — it would then constitute the 
standard t-J model for the ladder ]6). However, a closer 
look at the geometry of the CU2O5 coupled ladders reveals 
that even such a model would not be enough. In fact, one 
cannot neglect a Coulomb repulsion on-site element be- 
tween two holes on the same oxygen sites (oc U p 3> t p d 
0, where t p d is the Cu-O hybridization) but on two differ- 
ent oxygen 2p orbitals (x and y orbitals in Fig.QJ belong- 
ing to the two neighboring chains in different ladders Q. 
Thus, although the holes do not hop between the neighbor- 
ing CU2O5 ladders (a small inter-oxygen hopping can be 
neglected |8j) the t-J model for the single ladder should 
contain an extra interladder coupling which can be derived 
from this Coulomb interaction. 

Altogether this means that neglecting the kinetic and 
superexchange terms along the rung connecting the two 
chains and taking into account the interaction between the 
two neighboring chains is a priori as justified in the treat- 
ment of CU2O5 coupled ladders as skipping the interladder 
coupling and considering merely the t-J model for a sin- 
gle ladder. Naturally, the proper approach is to take into 
account both of these effects (under which the work is in 
progress) but in what follows we merely concentrate on the 
problem of two coupled chains connected by this interlad- 
der coupling, see Fig. Q] This will enable us to study the 
generic role of the interladder coupling. 

While the detailed derivation of the interladder cou- 
pling is somewhat lengthy and will be described elsewhere, 
a simple cartoon-picture shown in the left panel of Fig. [2] 
explains its basic idea. Here we see that the two Zhang- 
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Figure 2 Left panel: the artist's view of the interladder re- 
pulsion between two nearest neighbor Zhang-Rice singlets 
on two nearest neighbor chains. The Zhang-Rice state is 
depicted by a dotted ring. Large (small) arrows depict the 
hole spins for +1.0 (+0.25) charge. Right panel: the geome- 
try of the model described by the H+H+H' Hamiltonian, 
see text for further details. 



Rice singlets (£[ which are situated next to each other in the 
two neighboring chains (see dotted rings in Fig. |2]i share 
the same oxygen sites. More precisely, as an oxygen hole 
forming a Zhang-Rice singlet state is equally distributed 
among the four oxygen sites surrounding the central cop- 
per hole (9), there is a 0.25 hole charge in each of these 
oxygen orbital sharing the common oxygen sites. Since the 
two holes in two different orbitals repel as oc U p (where the 
Hund's exchange Jh is skipped for simplicity), the Zhang- 
Rice singlets repel as oc (1/4 x 1/4 + 1/4 x l/4)t/ p . Ob- 
viously, the Zhang-Rice states in each chain are not the 
proper orthogonal states (similarly as in the 2D case ||9l ) 
and one needs to orthogonalize them first. This, however, 
does not change the result qualitatively as it occurs that 
such a procedure yields that the two Zhang-Rice singlets 
repel as oc 0.14t7 p . Besides, finite Hund's exchange re- 
duces this interaction by a factor (1 — 5Jh/2U p ). 

Hence, we supplement the t-J model for a single 
chain, Eq. (G3, by the following interladder term which 
couples the two chains under consideration 



H = V 



(2) 



where V = 0.14:U p (1 — 5Jh/2U p ) and for realistic param- 
eters V ~ 0.5t. Here the bar signs are added over the 
operators which act in the Hilbert subspace of the neigh- 
boring chain. The detailed geometry of the full model for 
two coupled chains H = H + H + H' can be seen in 
the right panel of Fig. [2] — the solid lines show the bonds 
along which the kinetic and superexchange term of Eq. (Q} 
act, while the dashed lines show the bonds along which the 
interladder terms of Eq. (O are finite. 

3 Results In order to solve the t-J-V model, given 
by the Hamiltonian H, we first need to overcome the con- 
straint of no double occupancies. One of the approximate 



Copyright line will be provided by the publisher 



pss header will be provided by the publisher 



3 



ways to do it is to introduce the Gutzwiller factors g t = 
(2 — 2n) /(2 — n) and gj = 4/(2 — n) 2 , which renormalize 
the kinetic term (g t ) and the superexchange and interladder 
terms (<?j). Similar factors were used recently to investi- 
gate the superconducting flux phases in the cuprates iflOl . 
Then the whole Hamiltonian is defined in terms of unre- 
stricted fermions d\ a . 

Next, we decouple the fermion operators in the su- 
perexchange and interladder terms in a mean-field way: 

4a<k* d L' d i*' -> d L d ™( d L> d ™') + (4r d «r)4r'^' 

(a similar decoupling is used also for the bar operators). 
Then, we diagonalize the effective one-particle Hamilto- 
nian in k space assuming that the classical fields take their 
initial values as 



n — p for ij A € Z 
n + jtjP for ij A ^ Z 



(3) 



which defines the CDW order parameter p. As already 
mentioned, here A is the CDW period, cf. Fig. 1 in Ref. 
J3). In addition, (d\_i ^d^i^) are defined as in Eq. (f3]l 

but with i/X replaced by (i + 1)/A [(i + 2)/A for period 
A = 5] - such assumption minimizes the classical energy 
cost of the interladder repulsion V. Finally, the actual val- 
ues of the classical fields (and consequently the order pa- 
rameter p) are obtained by iterating the above procedure 
(i.e. calculated from the diagonalized Hamiltonian) until 
their initial and final values converge. 

As the main result of this self-consistent procedure, we 
obtained the increasing order parameter pas a function of 
the interladder interaction V for a realistic value of J = 
QAt |8| and for the three interesting hole dopings, see Fig. 
[3] First, we see that for finite values of V the CDW phase 
with the peculiar odd period is stable. Moreover, it is stable 
for realistic values of parameter V ~ 0.5t. Second, the 
CDW order with period A = 4 is less stable as the CDW 
order sets in only for V > 0.6i. 

One may ask, why the CDW phase with even period ap- 
pears to be less stable. Although this reproduces the exper- 
imental result quite well, it is a very peculiar phenomenon 
associated purely with the ID physics. In fact, it is the 
CDW phase with period A = 3 which is remarkably stable 
here due to the nesting of the band structure for this partic- 
ular doping, and indeed we find the exponential growth of 
p with increasing V > 0, special for the present ID model. 

4 Conclusions In summary, we showed that includ- 
ing the Coulomb repulsion between holes on the same oxy- 
gen sites but in different orbitals and belonging to two dif- 
ferent neighboring CU2O5 ladders, led to the effective re- 
pulsion between Zhang-Rice singlets which in turn led to 
the interladder repulsion V between nearest neighbor holes 
situated in the neighboring ladders. Next, we solved a t—J 
model supplemented by the repulsive term V in a simpli- 
fied geometry of two coupled chains. It occurred that the 
CDW state with odd period was stable in such a model for 
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Figure 3 The CDW order parameter p as a function of the 
interladder interaction V in the two coupled chains with 
J = QAt: for filling n = 2/3 (n h = 4/3) with period A = 
3 (solid line), for filling n = 3/4 (nh = 5/4) with period 
A = 4 (dashed line), and for filling n = 4/5 (rih = 6/5) 
with period A = 5 (dotted line). 



the realistic values of the interaction parameters V, while 
the CDW state with the even period (A = 4) was stable 
merely for slightly enhanced values of these parameters. 

The results presented in this paper could possibly ex- 
plain the onset of the novel odd-period-CDW state in the 
CU2O5 coupled ladder system in Sri4_ :r Ca. E Cu2404i with 
x = and x = 11. However, a detailed study of this phe- 
nomenon is currently under investigation since one needs 
to generalize the t—J—V model presented here to the cou- 
pled ladders and to investigate the stability of CDW states 
in this generalized model in order to verify this conjecture. 
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